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Ch ; Abstract 

A numerical study of the evolution of a massless scalar field in the background 
of rotating black holes is presented. First, solutions to the wave equation are 
obtained for slowly rotating black holes. In this approximation, the back- 
ground geometry is treated as a perturbed Schwarzschild spacetime with the 
angular momentum per unit mass playing the role of a perturbative param- 
■ eter. To first order in the angular momentum of the black hole, the scalar 

wave equation yields two coupled one-dimensional evolution equations for a 
function representing the scalar field in the Schwarzschild background and a 
qh second field that accounts for the rotation. Solutions to the wave equation are 

also obtained for rapidly rotating black holes. In this case, the wave equation 
does not admit complete separation of variables and yields a two-dimensional 
evolution equation. The study shows that, for rotating black holes, the late 
time dynamics of a massless scalar field exhibit the same power-law behavior 
as in the case of a Schwarzschild background independently of the angular 
momentum of the black hole. 
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I. INTRODUCTION 



Perhaps the simplest wave phenomenon in relativity is the propagation of linearized 
waves on a fixed, curved background. When observed at a fixed spatial point, the dy- 
namics of a wave propagating on spherically symmetric, time-independent, asymptotically 
flat, background geometries consists of three stages: During the first stage, the observed 
wave depends on the structure of the initial pulse and its reflection from the origin (Burst 
phase). This phase is followed by an exponentially decaying quasi- normal ringing of the 
black hole (Quasi- normal phase). In the last stage, the wave slowly dies off as a power-law 
tail (Tail phase). The last two phases are dictated by the interference of the part of the wave 
backscattered at the tail of the potential and the part reflected by the potential barrier. 

The tail phenomenon can be understood as due to the scattering of the wave off the 
effective curvature potential of the background geometry ^M. Tails have been mostly 
investigated on Schwarzschild backgrounds |3|-|5|], Reissner- Nordstrom black- holes ||, and 
for collapsing scalar fields [§,0. The study of tails has also implications in connection with 
the gravitational radiation emerging from inspiraling binary systems ||. 

For spherically symmetric systems, after separation of variables, the three-dimensional 
wave equation V M V M <I? = reduces, in a suitable radial coordinate x, to the one-dimensional 
wave equation: 

[-d? + dl-V(x)Mx,t) = 0, (1) 

where V(x) denotes the effective curvature potential. Recently, the late-time behavior of 
equation ([[]) has been the subject of detailed numerical and analytic investigations by Gund- 
lach et al. § and Ching et al. 

Gundlach et al. 0J^] studied both, the late-time behavior of massless fields around a fixed 
Schwarzschild geometry and the full nonlinear evolution of a minimally coupled scalar field. 
Their study shows a remarkable agreement between the numerically computed quasi-normal 
frequencies and power-law tails and their corresponding theoretical predictions. Using a null 
characteristic approach, they found tails not only at timelike infinity, but also at future null 
infinity and at the future horizon of the black hole. The time behavior of the tails is shown 
to be 

!^-{2l+p+x) a ^ timelike infinity, 

(t - x)~ il+pS) at future null infinity, (2) 
(t + x)~( 2l+p+lS> at future horizon. 

with P = 1 if there is an initial static field and P = 2 otherwise. 

Pursuing a different approach, Ching et al. showed that the late time tail phe- 
nomena are governed by the asymptotic structure of the background spacetime. Further- 
more, they found that, for non-static initial conditions, the simple power-law $ ~ £-22+0 
in Schwarzschild geometries, is not present for general potentials of the form V(x) ~ 
1(1 + l)/x 2 + x~ a In 2. Instead, the generic late-time behavior is given by $ ~ t- 2l+a Int. 

In spite of previous work, a complete understanding of the tail phenomenon for general 
systems is still missing. Questions such as which factors determine the magnitude and 
power-law time dependence have not been addressed for a general system. It is of particular 
interest, for instance, to show which role the dimensionality of the system plays on the 
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existence of power-law tails ||. When the point of view that tails arise from backscattering 
off an effective potential is adopted, it is not clear whether the tail behavior would remain 



unchanged, or be present at all, for intrinsically two or three dimensional systems JTC . 

A natural multi-dimensional generalization of the dynamics of scalar fields is obtained 
by considering the family of Kerr spacetimes. In contrast to the Schwarzschild background 
geometry, the wave equation V^V^ = in the background of a rotating black hole only 
admits separation of the azimuthal coordinate, so in principle a two-dimensional evolution 
problem has to be solved. It has been suggested ||, that tails similar to those present in 
Schwarzschild spacetimes should exist on Kerr backgrounds, a detailed analysis, however, 
has not been undertaken. 

This study considers the dynamics of scalar fields on the background of rotating black 
holes in two regimes: slowly and rapidly rotating holes. A considerable mathematical simpli- 
fication of the problem is achieved in the case of slowly rotating black holes since the angular 
momentum can be treated as a perturbative parameter. Thus, to first order in the angular 
momentum, the scalar wave equation yields two one- dimensional evolution equations. This 
formalism is described in Section [Dj. Results of the numerical calculations are presented in 
Section III], where power-law tails are discussed. The case of rapidly rotating black holes is 
discussed in Section [TV] 



II. SLOWLY ROTATING BLACK HOLES 



Using Boyer-Lindquist coordinates (t, r, 9, (f>) the wave equation for a massless scalar field 
$ reads 



(r 2 + a 2 ) 2 
A 



-a 2 sin 2 



_ _ AMar _ _ 
d n ® ^— d t ^ + 



1 



sin 2 9 A 



+d r (Ad r $) + — —do (sin 000$) = 
sin 6 



(3) 



where M is the mass of the black hole, a is its angular momentum per unit mass, and 
A = r 2 — 2Mr + a 2 . Equation (j3|) is equivalent to the special case of the Teukolsky equation 
fl^| for fields with vanishing spin weight. 

The Schwarzschild case, a = 0, allows separation of variables in terms of the scalar 
spherical harmonics YJ m (#, 0) without any requirements on a particular time behavior. The 
solution of (13T) can be written in the form 



$ = iyT(M)v|/(t,r*), 

r 



(4) 



where r* denotes the tortoise coordinate 

r* = r + 2M\n(r -2M) . (5) 
Substitution of (||) into (0), with a = 0, yields the one-dimensional wave equation 

\-d tt + d r * r , -V(r)]V(t,r*)=0, (6) 
with the potential V(r) defined by 
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It has been shown, in both numerical and analytic studies |1|J3|,§] , that at a fixed radius the 
solution of (HI) will fall off as t~( 2i + 3 ) for large t and non-static initial data. 

For a ^ 0, it is no longer possible to separate variables and arrive at an equation similar 
to @. While the azimuthal dependence of $ can still be described by e vm ^, a full separation 
of variables exists only if the time dependence is given by e lujt , as first demonstrated by Brill 
et al. In this case, both, the Teukolsky equation and equation (H), admit separable 



solutions of the form 



-iuit —im</> 



Sr(9;uj)R lm (r;w) 



(8) 



where gSj 11 are the spin- weight- s spheroidal harmonics. However, since the objective of 
the present study is the late-time dynamics, a decomposition based on u modes is not 
suitable, and one is forced, in principle, to solve a two-dimensional evolution equation for 
$. Nonetheless, the case of slowly rotating black holes circumvents the problem of solving 
a two-dimensional wave equation. 

The use of an azimuthal coordinate defined by asymptotic observers, as is the case with 
the Boyer-Lindquist <fi coordinate, introduces unphysical pathologies near the horizon even 
in the case a « M. A discussion of those pathologies and their precise manifestation in 
the slow rotation limit is given in the Appendix. Those coordinate induced problems are 
readily dealt with by adopting the Kerr azimuthal coordinate <fi given by 



d<f) = d<p + — dr . 

In the (t, r, 9, 0) coordinates, equation (|3]) is transformed into 

AMar 



(9) 



(r 2 + a 2 ) 2 
A 



— a sin 9 



sin 9 



dg(sin9d e <$>) + 



sin 2 fl 



d tt <5> 



A 

0, 



d r (Ad r $) 



(10) 



Under the assumption a « M, one can seek a solution of ( |I0D that possesses an angular 
dependence given by Y l m (9,4>) following a procedure similar to that used in deriving (f|). 
That is, one views the function ty(t,r*) containing a piece \Pq(£, r*) representing the scalar 
field in the Schwarzschild background and a second field \&i(t,r*) that takes into account 
the correction due to the rotation of the black hole to first order in a. Thus, the ansatz for 
the solution of ([10|) has the form 



(11) 



Before deriving the equations satisfied by and a further transformation of (10) to 
the coordinate system (t,r*,fi,<j>) is performed, where r* denotes the Schwarzschild tortoise 
coordinate @ and fi = cos 9. In these coordinates, substitution of (0) into the wave equation 
V^V^ = and collecting powers of a yields after separation of variables the equation 
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□o^ + a(no^i-po) + C(a 2 ) = 0, (12) 



where 

□o = -da + d r * r * - V(r) and (13a) 

-2Mr9 t f + r 2 a r .*o-(r-2M)* 
p Q = -2im — . (13b) 

The system of equations for \l/o an d ^i is obtained by requiring that the zero and first 
order terms of this expansion vanish independently, namely 

[-d u + d r * r * - V(r)\ ^o = 0, (14a) 

[-d tt + o r * r * + V[r)\ = 2™ (14b) 

with V(r) given by (0). As expected, the equation for does not contain terms depending 
on y&i, and the source term p in the equation for only depends on the solution \l/o of 
the zero order equation. The problem has then been reduced to solving the homogeneous 
equation (|14a|) for and using this solution as a source in the inhomogeneous equation 



(ED- 

Owing to the one-way membrane character of the horizon, the use of ingoing boundary 
conditions at the horizon arises most naturally for the solution of (|14a| ) . The asymptotic 
form of for r* — > — oo is then a wave with constant amplitude propagating to the horizon. 

As r — > 2M ( r* — > — oo ), the source term behaves as p — ► — d t ^o + d r *^/ . Thus, p - 
if "~ * ^o(t + r *)- On the other hand, at r* — > oo, p - due to the assumption that \l/ 
has compact initial data. 



III. POWER-LAW TAILS 

The numerical results presented here were computed with a second order staggered in 
time evolution scheme. For consistency, a null characteristic numerical integration was also 
used, and the results were in complete agreement. 

In general, both ty and are complex quantities. Owing to the purely imaginary 
coefficients of d t ty , d r *^ , and ^ in the source term in ( |14b| ), lm^ is coupled to Re\l/i, 



and Re\l/o to Im^i. Hence, without loss of generality, one may assume Re\l/o = Im^i = 0. 
Since the separation variable m plays the role of a scaling factor for the source term in ( 14b ), 
in the numerical calculations m = M = 1 is used. 

The initial data for \l/ consist of a bell-shaped pulse propagating outwards given by 

*!, m = ci [(c 2 (r* - 1) - r* n )(c 2 (r* - t) - r * ut )] 8 (15) 

for r in < r* < r out , and vanishing otherwise. The parameters have been chosen to yield a 
pulse centered at r* = 100; that is, c\ = 100, c 2 = 0.02 and r* nout = 1, 3. On the other hand, 
the initial conditions for are = dt^i = 0. Numerical experiments have shown, that 
the tail behavior is not affected by the choice of initial conditions for vE^. 

As previously discussed, the behavior of both the potential and the source term in ( |14b| ) 



makes it possible, in principle, to impose ingoing boundary conditions at the horizon and 
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outgoing conditions at infinity, i.e. lim r .^ TOO \l/ ,i = ^a,i(t ± r*). However, in a Cauchy 
evolution, such as the one here under consideration, boundary conditions are imposed at a 
Typically, the computational domains used covered —50 < r* < 2000 
0.1. Since r* = —50 yields r — 2M ~ 10~ 12 , imposing the left (ingoing into the 



finite distance |T5 
with Ar* 

black hole) boundary condition at a finite radius r* turns out to be a suitable approximation; 
both, the potential V{r) and the source po are negligible at r* = —50. A different situation is 
encountered at the right boundary. Even though the potential V(r) and the source po vanish 
as r* — > oo, at the right boundary (r* = 2000) of the computational domain, and especially 
for late times, the outgoing boundary condition is not a good approximation. However, this 
boundary is at a sufficiently large radial distance, so it allows enough dynamical evolution 
range to obtain the tail behavior before numerical boundary effects contaminate the solution. 

Figure [I] shows several snapshots of the evolution for / = 1, < t < 200. During the 
evolution, the initially outgoing zero order function ty is followed by a pulse in driven 
by the source p . In Fig. 1(a), the initial bell-shaped pulse for ^ centered at r* = 100 is 
shown. = po = 0. At t = 20, Fig. 1(b) depicts the radial dependence of the source p 
dominated by the initial shape of \I/o- As displayed in Fig. 1(c), after an evolution time of 
t = 100, a pronounced peak of po has developed due to the part of the initial pulse that was 
backscattered off the potential and is propagating inwards. In Fig. 1(d), at t — 200, as a 
consequence of the peak of p , a hump appears in the trailing part of 

Before presenting the numerical results, a heuristic analytical argument is given esti- 
mating the expected late time behavior. The argument is based on the result that the late 
time decay of waves propagating on curved spacetimes is dictated by the spatial asymptotics 
of the potential 0]. The argument does not include the contribution from the centrifugal 
barrier (/ = 0). Starting with the original wave equation ([!]) one obtains to first order in a 



-d tt + d r .r. - V{r)\ V + 



2imaM 



-2 d t + r d r * + 1 



2M\ 
r / 



$ = 0, 



(16) 



where 



$ = ^(t,r*)Y7™(#,0). 

Equation (|16|) consists of two operators: the operator for scalar fields on Schwarzschild back- 
ground (first term in squared brackets) and a second operator (first order in a) containing 
first order spatial and temporal derivatives. 

If one views the late-time tail behavior at a radial position r, as the result of the scattering 
by a potential at f >> r of a wave originated at r ~ r, the arrival time of the scattered 
wave is approximately given by t ~ (f — r ) + (f — r) « 2r At late times, the scalar field 
\& oc Vrif) w Vr(t), where Vx(r, t) = V(r) + V a (r, t) represents the "total" potential; that 
is, the Schwarzschild potential (^) with / = plus a correction V a due to the rotation of the 
black hole. This correction arises from the second term in squared brackets in ([!(]). Under 
the late-time and large scattering radius approximation, d t ^ ~ ^ /t and d r *^ ~ \l//r. Thus, 
V a (r,t) ~ f~ 3 + f~ 3 1" 1 + f~ 4 ~ t~ 3 . Therefore, \1/ oc V + V a ~ t~ 3 since V ~ t~ 3 . This shows 
that both, the Schwarzschild potential and the perturbative potential, contribute with the 
same power law behavior. The centrifugal barrier adds a t~ 2 1 factor to the tails H] . 

The numerical investigation of the late time behavior was performed at different distances 
from the black hole. Power-law behavior, with an exponent independent of the location of 
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observation, was detected for both \l/o and for different multipole indices I. As shown 
in Fig. and Fig. the numerically determined power-law exponents for the zero and first 
order solutions are in good agreement: For I = 1, the power-law exponents are —4.93 for 
\l/o an d —4.92 for NE^; the theoretically predicted value is 5 for this case. For I = 2, the 
exponents are —7.02 and —7.07, respectively, and the corresponding theoretical value is 7. 
These results lead to a power-law exponent for $ = - Y™ [\&o + a ^i] that is unchanged with 
respect to the behavior on a Schwarzschild background. 



IV. RAPIDLY ROTATING BLACK HOLES 

The next step is to consider the case of rotating black holes with arbitrary angular 
momentum per unit mass, i. e. < a < M. The objective here is to investigate whether 
the Schwarzschild tail behavior is modified for rapidly rotating black holes. The starting 
point is equation (|10"D written in coordinates (t,r* ,9, <ft), with the Kerr tortoise coordinate 
r* defined by 

dr* r 2 + a 2 



dr A 



^ = ^(t,r*,9)e im ^, (18) 



n _ (r 2 + a 2 ) 2 _ A ^ AimarM^^. 2 rA + % am (r 2 + a 2 )l 

-dtt^ - d r * r *y d ee ^> + d t ^> <9 r *^ 

a a a a 

COS ^ A TYl 2 A , . 

d e m + . ^ = 0, (19) 



Using the ansatz 



the wave equation (|10"D yields 



sin 6 1 cr sin ^ cr 

with 



cr = — r — r a 



4 ^2^2 i „^„2 



1 + cos 2 0) - 2a 2 rM sin 2 9 - a 4 cos 2 . (20) 

The wave equation (|T9"| ) reduces to the flat wave equation [— d u + d r * r *]^f = for r* — > ±oo, 
thus a simple formulation of radial outgoing boundary conditions is possible. 

Equation ( |T9"D was evolved on a rectangular spatial grid in r min < r < r max , < 9 < n/2. 
The numerical code was tested for convergence and stability. As in the one-dimensional case, 
the mass of the black hole was set to M = 1. The inner radial boundary r m j n was chosen so 
that the limit of the flat wave equation was reasonably approximated (cf. Section [TIT]). The 
data at 9 = 0,7r/2 were updated according to the behavior of the particular mode under 
consideration. The values of the field at the poles and the equatorial plane were determined 
by the value of m and the symmetry of the mode with respect to the equatorial plane: For 
m = 0, de^(t, r*, 0) = for all t, r*. For non-axisymmetric waves, given by m ^ 0, ^(0) = 
for all modes. At the equator, the condition for the modes symmetric to the equatorial plane 
is given by dgty(7r/2) = and the antisymmetric modes are characterized by \l/(7r/2) = 0. 

The initial data were given by Re\l/ = d t R.efy = and, using (|15D, Im^P = 



Pf a I yr 2 + a 2 , where a particular associated Legendre Polynomial P™ had to be chosen 
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according to the mode of interest. Because / does not appear in the equation, for fixed m, 
a mixing of modes belonging to different values of I will occur during the evolution if the 
initial data do not correspond to the lowest possible mode. 

In the following figures, the late-time behavior of at r* = 10 is displayed for different 
initial multipoles for an angular momentum parameter of a = 0.99. Fig. ^ shows the late 
time behavior for m — 1. Equidistant angular directions were chosen in the interval from 
9 = 7r/16 to 9 = 7r/2. The initial angular dependence corresponds to I — m — 1. The late- 
time behavior can be described by |^| ~ where the exponent is given by /x = 4.87 for 
all the displayed angles. In Fig. [5] the analogous situation is shown for m = 2 and an initial 
pulse corresponding to I = m = 2. In this case the power-law exponent for the observed 
angles is fi = 6.98. That is, the power-law exponents governing the late-time behavior for 
a = 0.99 do not exhibit a significant change when compared to the Schwarzschild case, in 
which the theoretical power-law tail exponents are given by fi = 5 for / = m = 1, and /i = 7 
for I = m = 2. In both cases the initial pulse used for the two-dimensional evolution is given 
by the lowest allowed mode for fixed m. 

The situation is different for the case I — 2, m = 0, depicted in Fig. || Here the initial 
pulse is not given by the lowest mode with m = that is symmetric with respect to the 
equatorial plane, corresponding to I = 0. Instead, the initial angular dependence was given 
by Im^ ~ P 2 ° ~ 3cos 2 # — 1. Here mixing of modes occurs and the late-time evolution is 
dictated by the lowest mode. The sink on the right-hand side is caused by the transition 
to the lowest mode. For large times the time dependence is given by \^\ ~ where 
2.88 < < 2.92 with < /i >= 2.91 when averaged over all observed angles. In contrast to 
this result, the corresponding Schwarzschild case exhibits no mixing of the modes, and the 
theoretical power-law tail exponent is given by /i = 7 for / = 2, m = 
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APPENDIX A: 

The horizon behavior of the wave equation written in Boyer-Lindquist coordinates is 
illustrated here, in the context of the slow rotation approximation. The manifestation of 
this behavior in the solution of the initial value problem is analysed in the limit r — > 2M. 

The small a expansion of equation (|3|) using the ansatz 

$ = - r Y t m (9, <P) [* (t, r*) + a ^(t, r*)} , (Al) 
yields the system of equations 

[-d tt + d r * r * - V(r)] ^ = (A2a) 

Am M 

[-du + <W - V(r)} d t ^ . (A2b) 



S 



Upon examining equation ( |A2b|) a problematic feature emerges. The source appearing 



on the right hand side does not vanish in the limit r* — > — oo. More precisely, the asymptotic 
form of the equation at the horizon is (— d u + <9 r . r *)^i = jffz dt^o- Using null coordinates 
v = t + r* and u = t — r*, and assuming a purely ingoing solution ^o(v), one obtains, 

^(u,v) = b(u)Mv) + fi(v) + h{u) , (A3) 

where 

b{u) = -MP U > (A4) 

and fx, fn are solutions of the homogeneous equation. 

The presence of this growing mode in the system may at first appear puzzling, since the 
governing equations @ and fllOD are related simply by an angular coordinate redefinition. 
Yet the singular nature of the transformation (|D near the horizon is the source of the growing 
mode. A direct comparison of systems (|A2|) and (|14]) reveals that the two sets of equations 
are equivalent, with given by 

i?x{t,r*) = #i(*,r*) + c(r)* (t,r*) , (A5) 

where 

, , ixm im , , r s ... 

c ' r » = 2M-2M ln <7^)- < A6 > 



The solution of the equation with the correct source given by ( |13b|) is denoted by 



r 



The coefficient c(r) is singular at the horizon, hence bounded solutions ^/q, give rise 
to an unbound combination for The linear combination ( |A5| ) has indeed the horizon 
singularity demonstrated by (|A3|) . This singularity follows from the fact that 



lim c(r) = lim b(u) pa lim ln(r — 2M) . (A7) 

An intuitive explanation of this pathology lies in the twisting of azimuthal directions defined 
at infinity, as seen by an infalling observer near the horizon. The wave evolution is well 
behaved in local frame, yet it appears singular in a coordinate system that winds itself 
infinitely many times around the black hole. The singularity of the evolution in the Boyer- 
Lindquist coordinates is introduced entirely through this singular (at the horizon) coordinate 
transformation. 
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FIGURES 
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FIG. 1. Evolution of'tfo, Po, and *i for for I = 1, < t < 200. (a) At t = 0, *i = p = 0. 
(b) At t = 20, the source po is dominated by the initial shape of ^o- (c) At t = 100, a peak in po 
has developed due to the part of the initial pulse that was backscattered off the potential and is 
propagating inwards, (d) At t = 200, as a consequence of the peak of po in (c), a hump appears in 
the trailing part of 
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FIG. 2. Log-log plots of and *i for I = 1 at r* 
—4.93 for \I>o and —4.92 for The wiggles that can be seen for t < 200 are remainders of the 
quasi-normal ringing of the black hole. 
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FIG. 3. Log-log plots of and for I = 2 at r* = 10. The power-law exponents are given 
by -7.02 and -7.07, respectively. 
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FIG. 4. Log-log plots of |^| for m = 1, a = 0.99 at r* = 10 for different equidistant angular di- 
rections from 6 = 7r/16 to 6 = tt/2. The initial angular dependence was given by Im* ~ P/ ~ sin 9. 
For large times the time dependence is given by |^| ~ t~^, where /U = 4.87 ± 0.005 for all observed 
angles. The theoretical power-law tail exponent in the Schwarzschild case is given by fj, = 5. 
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FIG. 5. Same as in Fig. 4, but for m = 2 and the initial angular dependence given by 
Im^ ~ P| ~ For large times the power-law exponent is \i = 6.98 ± 0.004 for all observed 

angles. The theoretical value in the Schwarzschild case is given by /j, = 7. 




FIG. 6. Log- log plots of j^l = ImVl/ for m = 0, a = 0.99 at r* = 10 for different angular 
directions in equidistant steps from 9 = to 9 = ir/2. The initial angular dependence was given 
by Im* ~ P 2 ° ~ 3 cos 9 — 1. Mixing of modes: For large times the time dependence is given 
by |^| ~ t~^, where for Im* the exponent is 2.88 < /(/ < 2.92 with < fi >= 2.91 over the 
observed angles. Re^ = at all times, because the two equations are decoupled for m = 0, and 
Re^ = 4rRe$> = initially. 



14 



